Quasiequilibrium sequences of black-hole— neutron-star binaries in general relativity 
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We construct quasiequilibrium sequences of black hole-neutron star binaries for arbitrary mass 
ratios by solving the constraint equations of general relativity in the conformal thin-sandwich decom- 
position. We model the neutron star as a stationary polytrope satisfying the relativistic equations 
of hydrodynamics, and account for the black hole by imposing equilibrium boundary conditions 
on the surface of an excised sphere (the apparent horizon). In this paper we focus on irrotational 
configurations, meaning that both the neutron star and the black hole are approximately nonspin- 
ning in an inertial frame. We present results for a binary with polytropic index n = 1, mass ratio 



5 and neutron star compaction M^dm,o/^o = 0.0879, where is the irreducible 



mass of the black hole, Afg the neutron star baryon rest-mass, and AfADM.o ^nd _Ro the neutron 
star Arnowitt-Deser-Misner mass and areal radius in isolation, respectively. Our models represent 
valid solutions to Einstein's constraint equations and may therefore be employed as initial data for 
dynamical simulations of black hole-neutron star binaries. 

PACS numbers: 04.30.Db, 04.25. Dm, 04.40.Dg 



Coalescing black hole-neutron star (hereafter BHNS) 
binaries are among the most promisingsources of gravi- 
tational waves for laser interferometers [ll Qi S U ■ BHNS 
mergers may reveal a wealth of astrophysical informa- 
tion (see e.g. 0), and, along with mergers of binary 
neutron stars, are also considered primary candidates 
for central engines of short-duration gamma-ray bursts 
(SGRBs) 11013. Recent observations of several SGRBs 
localized by the Swift and HETE-2 satellites in regions 
with low star formation strongly suggest that a compact 
binary merger scenario for SGRBs is favored over models 
involving the collapse of massive stars (see, e.g., and 
references cited therein). 

Significant effort has gone into the study of binary 
neutron stars and binary black holes, which are also 
promising sources of gravitational radiation. Fully rel- 
ativistic simulations of BHNS binaries have received 
far less attention. Most BHNS calculations to date, 
including quasiequilibrium (QE) calculations 0, ITll 
11^ [1^ Im IT^ \wL [tR IT^ and dynamical treatments 
[T9ll20l l2l L l22 l l23 ll24l25l |. employ Newtonian gravitation 
in either some or all aspects of their formulation. We have 
recently launched a new effort to study BHNS binaries 
in a fully relativistic framework (see also [26L l27l|). first 
by constructing QE models I23 and then by employ- 
ing them as initial data in dynamical simulations ^ j30] . 
So far we have focused on binaries for which the black 
hole mass is much greater than the neutron star mass. 
For binaries with such extreme mass ratios the rotation 
axis can be taken to pass through the center of the black 
hole, and the tidal effects of the neutron star on the black 
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hole may be ignored. These app roximations simplify the 
problem considerably (see |28l|). However, they break 
down for binaries containing comparable mass compan- 
ions. Such systems are more suitable as SGRB candi- 
dates, because the tidal disruption of the neutron star by 
the black hole will occur near or outside the innermost 
stable circular orbit. This disruption may be necessary 
to create a gaseous accretion disk around the black hole 
capable of generating a SGRB 0- Gravitational waves 
from BHNS binaries of comparable mass are detectable 
by ground based laser interferometers like LIGO (Laser 
Interferometric Gravitational- wave Observatory), while 
waves from systems with extreme mass ratios are much 
lower in frequency and require space-borne interferome- 
ters like LISA (Laser Interferometer Space Antenna). 

In this paper we describe the construction of QE se- 
quences of BHNS binaries with companions of compa- 
rable mass. We construct such binaries by solving the 
constraint equations of general relativity together with 
the relativistic equations of hydrodynamic equilibrium 
in a stationary spacetime assuming the presence of an 
approximate helical Killing vector (see, e.g., the recent 
reviews |31|, as weU as Sec. II of ,3^). Throughout 
this paper we adopt geometric units with G = c = 1, 
where G denotes the gravitational constant and c the 
speed of light. Latin and Greek indices denote purely 
spatial and spacetime components, respectively. 

The line element in 3-f 1 form can be written as 



= -a^df + 7y [dx' + (3'dt){dx^ + /3^dt), (1) 



where a is the lapse function, the shift vector, 7^ the 
spatial metric, and g^i, the spacetime metric. Einstein's 
equations can then be split into constraint and evolution 
equations for the spatial metric "fij. To decompose the 
constraint equations we introduce a conformal rescaling 
7ij — '^'^lij J where "0 is the conformal factor and the 



2 



spatial background metric. The Hamiltonian constraint 
then reduces to 

V^V = -27r^V + I^JR + l^V-'^' - l^J-'A,,A'\ (2) 
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Here V;, Rij, and R — Rij denote the covariant 
derivative, the Ricci tensor, and the scalar curvature as- 
sociated with 7,y. We also decompose the extrinsic cur- 
vature K''^ into its trace {K) and traceless {A^-') parts. 

In the conformal thin-sandwich decomposition we ex- 
press the traceless part of the extrinsic curvature in 
terms of the time derivative of the background metric, 
Uij = dt"fij, and gradients of the shift. For the construc- 
tion of equilibrium data it is reasonable to assume Uij = 
in a corotating coordinate system, which yields 

A^i = p (v*/?^" + V^"/3' - If^Vkf^") . (3) 

Inserting this into the momentum constraint we then ob- 
tain 

3 

= WTTaij'^f + 2i*^ Vj(aV"®) + -af^VjK. (4) 
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It is also reasonable to assume dtK = 0, which, from the 
evolution equation for the extrinsic curvature, yields 

V^a = 4:TTai}^{p + S) + ai}-^A,jA'^ - 

2fiV,aVj \aij*K^ + tp*P'V^K. (5) 

Before we can solve the above set of gravitational field 
equations for ip^ /?* and a, we still need to specify the 
spatial background metric 7ij and the trace of the extrin- 
sic curvature K. We choose this background geometry 
to describe the Schwarzschild metric expressed in Kerr- 
Schild coordinates. Specifically, we choose jij = rjij + 

2MBiihlj/rBH and K = 2d|HMBH(l + 37\fBH/rBH)/r-|H- 
Here 77^ is the flat spatial metric, Mbh is the "bare" mass 

of the black hole, tbh = (^|h + ^bh + ^bh)^^^ is the 
coordinate distance from the black hole center, k = P = 
JTgjj/rBH is the radial vector pointing away from the 
black hole center, and Abh = (1 + 2Mbh/''bh)^^^^ is the 
lapse function of the Schwarzschild metric in Kerr-Schild 
coordinates. The matter terms on the right-hand side of 
Eqs. l(2Jl, Q, and ||SJ| are the projections p = n^n^T^^''' ^ 
f = -7>.r^^ = ir^jj^T'^", and S = -f'^ of the 
stress-energy tensor T^,y, where is the unit vector nor- 
mal to the spatial hypersurface. Assuming an ideal fluid 
we have T^,, = (po + Pi + P)u^Ui, -\- Pg^,^, where is 
the fluid 4-velocity, po the baryon rest-mass density, pi 
the internal energy density, and P the pressure. 

The elliptic Eqs. ©, (gj), and © require boundary 
conditions, both at spatial inflnity and on the surface of 
an excised sphere within the black hole interior. At spa- 
tial infinity, where the metric becomes asymptotically flat 



in an inertial frame, we impose the exact boundary condi- 
tions, and on the excision surface (apparent horizon) we 
impose the black hole equilibrium boundary conditions 
suggested in '34]. To construct approximately nonspin- 
ning black holes we set the shift according to Eqs. (39) 
and (50) in |34j with Q,. = ilo in their notation. In 
the language of fs^l this assignment corresponds to the 
"leading-order approximation" , and we plan to improve 
this approximation as outlined there. 

In addition to the field Eqs. 101, Q, and (0) we have 
to solve the equations of relativistic hydrodynamics. For 
stationary configuration the relativistic Euler equation 
can be integrated once to yield 

haj/jo — constant, (6) 

where h = {po + pi + P)/ po is the fluid specific enthalpy, 
and 7 and 70 are Lorentz factors between the fluid, the 
rotating frame, and the inertial frame (see Sec. II. C. of 
''29'| for the definitions). For irrotational fluids the fluid 
velocity can be expressed in terms of the gradient of a 
velocity potential The equation of continuity then 
becomes 

{po/h)WV^^ + (V^vE')V^(po/M = 0, (7) 

where is the covariant derivative associated with g^^. 
We solve these equations for a polytropic equation of 
state P = kpq , where F = 1 -I- 1/n denotes the adiabatic 
index, n is the polytropic index and k is a constant. Here, 
we focus on the case n = 1 (i.e., F = 2). 

We determine the orbital angular velocity by requiring 
that the derivative in the X direction of the enthalpy field 
at the center of the neutron star be zero, which implies 
that the total force balances at the center of the neutron 
star [23 . We confirm that the angular velocity obtained 
by this method agrees with that obtained by requiring 
the enthalpy at two points on the neutron star's surface 
be equal to within one part in 10~^ "2^. 

We locate the axis of rotation by requiring that the to- 
tal linear momentum P* ^ 'g^ §00 K^^dSj vanish To 
do so, we first align the axis of rotation with the Z axis 
and place the X axis to be the perpendicular line to the Z 
axis that passes through the black hole center. Given the 
equatorial symmetry in the problem, the Z component of 
the momentum vanishes automatically. With the orbital 
angular velocity held fixed, we drive the Y component 
of the linear momentum toward zero by adjusting the X 
coordinate of each companion, keeping their separation 
in the X direction unchanged. To determine their Y co- 
ordinates, we require that the X component of the linear 
momentum be zero, but only adjust the Y coordinate of 
the neutron star to achieve this. Thus, we fix the black 
hole's center to remain on the X axis at F = 0. 

Our numerical code uses the spectral method Lorene 
library routines developed by the Meudon relativity 
group 37]. The computational grid is divided into 10 (8) 
domains for a black hole (neutron star) and its exterior, 
and each domain is covered by TV,- xNgX = 33 x 25 x 24 
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TABLE I: Physical parameters for a binary sequence with 
mass ratio M^^/M^^ = 5 and neutron star compaction 
M^DM o/^o = 0.0879 (where Ro is the areal radius of the iso- 
lated neutron star). The baryon rest mass, the ADM mass, 
and the isotropic coordinate radius of the neutron star in iso- 
lation are M^^ = 0.1, M^dm.o = 0.0956, and fo = 0.990 
(k = 1). We list the binding energy i?b, total angular mo- 
mentum J, orbital angular velocity fi, maximum density pa- 
rameter (/max, mass-shedding indicator Xmin, and fractional 
difference 5M between the ADM mass and the Komar mass. 



d/Mo E^/Mg JjMl mio 5M 

20.46 -7.18(-3) 0.679 1.10(-2) 5.83(-2) 0.892 7.95(-4) 

18.41 -8.11(-3) 0.648 1.28(-2) 5.84(-2) 0.913 3.12(-3) 

16.36 -9.22(-3) 0.615 1.52(-2) 5.86(-2) 0.953 6.41(-3) 

14.32 -1.03(-2) 0.586 1.86(-2) 5.87(-2) 0.895 1.00(-2) 

12.29 -1.13(-2) 0.556 2.36(-2) 5.75(-2) 0.843 1.43(-2) 

10.26 -1.34(-2) 0.520 3.07(-2) 5.61(-2) 0.845 1.89(-2) 

9.243 -1.47(-2) 0.502 3.59(-2) 5.56(-2) 0.791 2.07(-2) 

8.741 -1.53(-2) 0.494 3.91(-2) 5.48(-2) 0.710 2.19(-2) 

8.439 -1.58(-2) 0.488 4.11(-2) 5.42(-2) 0.588 2.27(-2) 
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FIG. 1: Contours of the lapse a in the equatorial plane for the 
innermost configuration of the sequence listed in Table Q The 
thick circle on the left denotes the excision surface (apparent 
horizon) of the black hole, while that on the right denotes 
the surface of the neutron star. The cross "x" indicates the 
position of the rotation axis. 



collocation points except for the closest two where 9 (8) 
domains are covered by A^r x A^e x A^^ = 25 x 21 x 20 
points. 

In the following we focus on results for an inspiral 
sequence of constant irreducible black hole mass 
and neutron star baryon rest mass 

M]^S Table n we 
list results for a binary of mass ratio Mi^H/M^S = 5 
and a neutron star mass Afg^ =0.1, where the bar de- 
notes non-dimensional polytropic units M — k~^^'^M. 
At infinite separation, this neutron star has a compaction 
-^^ADM o/^o = 0.0879, where Rq is the areal radius of the 
spherical star in isolation. In Fig.^we also show contours 
of the lapse a for the innermost configuration of this se- 
quence. The value of the lapse everywhere on the black 
hole apparent horizon is set to be its Kerr-Schild value 
there, 2~^/^ = 0.7071; its value at the center of the neu- 
tron star is 0.7275. Note that the center of the neutron 
star, as defined by the maximum of the enthalpy, does 
not coincide with the minimum point of the lapse inside 
the star. These configurations are the first fully rela- 
tivistic QE models of BHNS binaries that do not assume 
an extreme mass ratio and employ equilibrium boundary 
conditions to model the black hole. 

In the table we list the fractional binding energy 
E\j/Mi) = Madm/-^^o ~ Ij total angular momentum 
J, orbital angular velocity fi, maximum of the density 
parameter (/max = (-P/po)max, minimum of the mass- 
shedding indicator x = (9(lnft,)/9r)cq/(i9(ln/i)/9r)poic 
[2^ im , and fractional difference between the Arnowitt- 
Deser-Misner (ADM) mass and Komar mass 5M = |1 — 
AfKom/-MADM|- Quantities are tabulated as functions of 
the coordinate separation between the center of the black 
hole and the point of maximum baryon rest-mass density 



in the neutron star. Here, A/q = + M^y)m ^^e 

ADM mass of the binary system at infinite orbital separa- 
tion, i.e., the sum of the irreducible mass of the isolated 
black hole and the ADM mass of an isolated neutron 
star with the same baryon rest-mass. For an isolated 
Schwarzschild black hole, the ADM mass is the same as 
the irreducible mass. 

We compared our values for the angular momentum 
with those from third-order post-Newtonian (3PN) ap- 
proximations [33 and found agreement to within about 
5% for close binaries, and better agreement for larger 
separations. In agreement with the 3PN results we do 
not see any indication of a turning point in the angular 
momentum, meaning that the tabulated sequence does 
not exhibit an innermost stable circular orbit, hence the 
binary orbits are all stable. 

The quantity Xmin is defined by the minimum of the 
indicator x which compares the gradient of In h at the 
pole with that on the equator. For spherical stars at 
infinite separation we have Xmin = 1, while Xmin = in- 
dicates the formation of a cusp and hence tidal breakup. 
As we have discussed in , spectral methods no longer 
converge in the presence of discontinuities, so that our 
sequence terminates before reaching Xmin ~ 0. However, 
extrapolating from the last three data points we estimate 
that the star will be tidally disrupted when fiAfo ~ 0.046. 
This value agrees with those estimated via the approx- 
imate relativistic expansion of llSll (f^Afn ~ 0.043) and 
the purely Newtonian models of [13 (ilAfo ~ 0.046). 

Equality between the ADM mass and Komar mass is 
equivalent to satisfying a relativistic virial theorem and 
indicates that the system is stationary (cf. js^)- In our 
calculation we do not impose this equality to construct 
the sequence, but instead evaluate mass difference as a 
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diagnostic. For decreasing separation the fractional dif- 
ference between the two masses increases to over 2% for 
our innermost configuration. This clearly indicates that 
our closest models are not in perfect equilibrium. The re- 
sulting small, but finite, systematic mass difference has a 
large effect on the binding energy, which is also computed 
as the small difference between much larger masses. 

We speculate that the differences between the ADM 
and Komar masses could be caused by our choice of the 
background geometry. Our choice of a Kerr-Schild back- 
ground metric is motivated by our requirement that it 
correctly reduce to the exact solution for a spinning black 
hole in the limit of large separation (although here we 
only treat nonspinning holes). Also, in our coordinates, 
the lapse remains positive on the horizon ( "horizon pen- 
etration"), which is necessary when computing Aij from 
Eq. (PJ. In we compared with a flat background 
and found that the choice of the background has a small 
but non- vanishing effect on the physical properties of the 
resulting binary configurations (see also the discussion of 
non-maximally sliced black hole binaries in [3^1. moti- 
vating our speculation that our choice here may result in 
the small but systematic deviation from perfect equilib- 
rium. Our configurations are solutions to the constraint 
equations, and are hence adequate initial data for dy- 
namical simulations of BHNS binaries (our main moti- 
vation here). We are currently experimenting with other 
background solutions to find better approximations to 
quasiequilibrium. However, this discrepancy in mass also 
heightens interest in the recent "wave-less" formulation 
of the initial value problem that is based on the equality 
of the ADM and Komar masses and avoids the need to 
choose a background geometry altogether ,4Q„ ^] . 

Finally, we turn our attention to the validity of QE 
configurations in circular orbit as initial data for dynam- 
ical simulations. The assumption of circular orbits and 
an associated helical Killing vector for relativistic bina- 
ries is an approximation, since the emission of gravita- 
tional waves leads to orbital decay. This approximation 
breaks down at a certain binary separation when the in- 
spiral can no longer be ignored. We can quantify the 
departure from true QE by comparing the time-scale of 



the orbital period with that of the orbital decay driven 
by the emission of gravitational waves. To lowest order 
we can estimating the ratio between these two timescales 
with the help of the quadrupole formula for Newtonian 
point-masses, which yields 




Here 1^ = M^^M^^ /{M^^+M^^)'^ and d^in denotes the 
closest binary separation we computed in this paper (the 
value of the last line in Table QJ. It is reasonable to ap- 
proximate the binary orbit as circular as long as the ratio 
^orb /tow is si gni ficantly smaller than unity (compare the 
discussion in |42]). For d = c?min, we have torb ~ 0.2tGW, 
and for larger separations the ratio torb /tow isXls off with 
^-5/2 Pqj. ^jjggg separations it is therefore reasonable to 
neglect the inspiral and construct binaries in circular or- 
bits and in the presence of a helical Killing vector. 

In summary, we compute sequences of BHNS bina- 
ries with comparable mass companions. We solve the 
constraint equations of general relativity in the confor- 
mal thin-sandwich decomposition, subject to equilibrium 
black hole boundary conditions, together with the rel- 
ativistic equations for hydrodynamic equilibrium in a 
stationary spacetime. We construct irrotational bina- 
ries, adopt a polytropic equation of state for the neu- 
tron star, and choose the background geometry to be a 
Schwarzschild black hole expressed in Kerr-Schild coor- 
dinates. As an example, we present results for a binary 
of mass ratio M^^/M^^ = 5 and neutron star of com- 
paction Mf^y^ J/Ro = 0.0879. To the best of our knowl- 
edge, these are the first models of quasiequilibrium, cir- 
cular orbit, relativistic BHNS binaries with companions 
of comparable mass. 

J.A.F. is supported by NSF Grant AST-0401533. 
This paper was supported in part by NSF Grants 
PHY-0205155 and PHY-0345151, and NASA Grant 
NNG04GK54G, to University of Ilhnois at Urbana- 
Champaign, and NSF Grant PHY-0456917 to Bowdoin 
College. 



[1] G. Gonzalez (LIGO Scientific Collaboration), Class. 

Quantum Grav. 21, S691 (2004). 
[2] M. Hewitson, P. Aufmuth, C. Aulbert, S. Babak, and 

R. Balasubramanian, Class. Quantum Grav. 20, S581 

(2003). 

[3] M. Ando (TAMA Collaboration), Class. Quantum Grav. 

19, 1409 (2002). 
[4] F. Acernese (VIRGO Collaboration), Class. Quantum 

Grav. 21, S709 (2004). 
[5] M. Vallisneri, Phys. Rev. Lett. 84, 3519 (2000). 
[6] M. Shibata and K. Taniguchi, Phys. Rev. D 73, 064027 

(2006). 

[7] J. A. Faber, T. W. Baumgarte, S. L. Shapiro, and K. 



Taniguchi, Astrophys. J. 641, L93 (2006). 
[8] D. J. Price and S. Rosswog, accepted to Science (2006), 

astro-ph/0603845 
[9] E. Berger, appearing in Proceedings of the 16th Annual 

October Astrophysics Conference in Maryland, "Gamma 

Ray Bursts in the Swift Era", eds. S. Holt, N. Gehrels, 

and J. Nousek, astro-ph/0602004 
[10] S. Chandrasekhar, Ellipsoidal Figures of Equilibrium 

(Yale University Press, New Heaven, CT, 1969). 
[11] L. G. Fishbone, Astrophys. J. 185, 43 (1973). 
[12] D. Lai, F. A. Rasio, and S. L. Shapiro, Astrophys. J. 

Suppl. Ser. 88, 205 (1993). 
[13] D. Lai and A. G. Wiseman, Phys. Rev. D 54, 3958 (1996). 



5 



[14] K. Taniguchi and T. Nakamura, Prog. Theor. Phys. 96, 
693 (1996). 

[15] M. Shibata, Prog. Theor. Phys. 96, 917 (1996). 
[16] K. Uryu and Y. Eriguchi, Mon. Not. R. Astron. Soc. 303, 
329 (1999). 

[17] P. Wiggins and D. Lai, Astrophys. J. 532, 530 (2000). 
[18] M. Ishii, M. Shibata, and Y. Mino, Phys. Rev. D 71, 

044017 (2005). 
[19] B. Mashhoon, Astrophys. J. 197, 705 (1975). 
[20] B. Carter and J.-P. Luminet, Astron. Astrophys. 121, 97 

(1983); Mon. Not. R. Astron. Soc. 212, 23 (1985). 
[21] J.-A. Marck, Proc. R. Soc. Lond. A 385, 431 (1983). 
[22] W. H. Lee and W. Kluzniak, Astrophys. J. 526, 178 

(1999); Mon. Not. R. Astron. Soc. 308, 780 (1999). 
[23] W. H. Lee, Mon. Not. R. Astron. Soc. 318, 606 (2000); 

328, 583 (2001). 
[24] S. Rosswog, R. Speith, and G. A. Wynn, Mon. Not. R. 

Astron. Soc. 351, 1121 (2004). 
[25] S. Kobayashi, P. Laguna, E. S. Phinney, and P. Meszaros, 

Astrophys. J. 615, 855 (2004). 
[26] M. Miller, gr-qc/0106017 (unpublished). 
[27] C. F. Sopuerta, U. Sperhake, and P. Laguna, to appear 

in Class. Quantum Grav. (2006), gr-qc/0605018 
[28] T. W. Baumgarte, M. L. Skoge, and S. L. Shapiro, Phys. 

Rev. D 70, 064040 (2004). 
[29] K. Taniguchi, T. W. Baumgarte, J. A. Faber, and S. L. 

Shapiro, Phys. Rev. D 72, 044008 (2005). 



[30] J. A. Faber, T. W. Baumgarte, S. L. Shapiro, K. 
Taniguchi, and F. A. Rasio, Phys. Rev. D 73, 024012 
(2006). 

[31] G. B. Cook, Living Rev. Relativity 3, 5 (2000). 

[32] T. W. Baumgarte and S. L. Shapiro, Phys. Rep. 376, 41 

(2003) . 

[33] E. Gourgoulhon, P. Grandclement, K. Taniguchi, J.-A. 
Marck, and S. Bonazzola, Phys. Rev. D 63, 064029 
(2001). 

[34] G. B. Cook and H. P. Pfeiffer, Phys. Rev. D 70, 104016 

(2004) . 

[35] M. Caudill, G. B. Cook, J. D. Grigsby, and H. P. Pfeiffer, 
submitted to Phys. Rev. D, gr-qc/0605053 

[36] J. M. Bowen and J. W. York, Jr., Phys. Rev. D 21, 2047 
(1980). 

[37] LORENE web page (http:/ /www. lorene.obspm.fr /J. 

[38] L. Blanchet, Phys. Rev. D 65, 124009 (2002). 

[39] E. Gourgoulhon, P. Grandclement, and S. Bonazzola, 

Phys. Rev. D 65, 044020 (2002). 
[40] M. Shibata, K. Uryu, and J. L. Friedman, Phys. Rev. D 

70, 044044 (2004).' 
[41] K. Uryri, F. Limousin, J. L. Friedman, E. Gourgoulhon, 

and M. Shibata, submitted, gr-qc/0511136 
[42] M. D. Duez, T. W. Baumgarte, S. L. Shapiro, M. Shibata, 

and K. Uryu, Phys. Rev. D 65, 024016 (2002). 



